Optimal Light Beams and Mirror Shapes for Future LIGO Interferometers 
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We report the results ol a recent search for the lowest value of thermal noise that can be achieved 
in LIGO by changing the shape of mirrors, while fixing the mirror radius and maintaining a low 
diffractional loss. The result of this minimization is a beam with thermal noise a factor of 2.32 (in 
power) lower than previously considered Mesa Beams and a factor of 5.45 (in power) lower than the 
Gaussian beams employed in the current baseline design. Mirrors that confine these beams have been 
found to be roughly conical in shape, with an average slope approximately equal to the mirror radius 
divided by arm length, and with mild corrections varying at the Fresnel scale. Such a mirror system, 
if built, would impact the sensitivity of LIGO, increasing the event rate of observing gravitational 
waves in the frequency range of maximum sensitivity roughly by a factor of three compared to an 
Advanced LIGO using Mesa beams (assuming all other noises remain unchanged). We discuss the 
resulting beam and mirror properties and study requirements on mirror tilt, displacement and figure 
error, in order for this beam to be used in LIGO detectors. 
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First-generation laser interferometer gravitational- 
wave detectors, such as initial LIGO [l[, VIRGO Q, 
GEO 600 @ and TAMA 300 0, have either reached or 
approached their respective design sensitivities, and have 
taken a first round of coordinated scientific data. While 
a detection in the near future is possible with these first 
data, or with upcoming data from Enhanced LIGO, a 
moderate upgrade of LIGO detectors, the construction of 
second generation detectors, e.g., Advanced LIGO, which 
are at least 100 times more sensitive than first-generation 
detectors (in power), has already started. 

Advanced LIGO reaches its maximum sensitivity in 
the frequency range 50-300 Hz. In this band, the inter- 
nal thermal noise, i.e., thermal fluctuation in the mean 
location of the mirror's surface of reflection, relative to 
the mirror's center of mass, is the dominant noise source. 
Depending on location, internal thermal noise can be 
divided into coating thermal noise and substrate ther- 
mal noise; while from the physical origin, noise that arise 
from internal friction (viscosity) is called Brownian noise, 
while noise that corresponds to thermal damping is called 
thcrmoelastic noise. The dominant component of in- 
ternal thermal noise in fused silica mirrors, the leading 
choice for Advanced LIGO, is the noise contributed by 
the mirror coating 

Lowering internal thermal noise will not only directly 
increase LIGO's event rate and thus our chances of see- 
ing gravitational waves [29], but it may also help bring 
Advanced LIGO sensitivity beyond the Standard Quan- 
tum Limit. Thus, for the first time, LIGO can study 
uantum effects as experienced by 40-kilogram objects 
(Initial LIGO will not benefit from lowering inter- 
nal thermal noise, because in this frequency band it will 
be limited by shot noise and suspension thermal noise.) 
Other advanced ground-based detectors as well as Quan- 
tum Non-Demolition experiments 0, Q will also benefit 
from implementing the research described here when bat- 
tling mirror internal thermal noise. 



Advanced LIGO's present design (baseline design) uses 
arm cavities with Gaussian light beams supported by 
spherical mirrors — with waist-size of Gaussian beams 
chosen to maintain a ~ 1 ppm per bounce diffrac- 
tional loss. As proposed by Vinet Q and Thorne 
and O'Shaughnessy [9(, wider beams supported by non- 
spherical mirrors can average better over thermal fluctu- 
ations throughout the mirror surface, and therefore give 
less thermal noise than conventional cavities with spher- 
ical mirrors and Gaussian beams. As a straightforward 
example, Thorne and O'Shaughnessy proposed the so- 
called Mesa beam, which can lower coating noise by a 
factor of 2.35 0]. 

In this work, we aim to search for the beam/mirror 
configuration with the minimum possible thermal noise. 
In order to do so, a systematic search over all possible 
freely propagating beams will be performed. More pre- 
cisely, to ensure the completeness of the search, we ex- 
pand the fundamental optical mode of the cavity using 
a Gauss-Laguerre basis (with real-valued coefficients) at 
the center of the cavity, while defining mirror shapes by 
the beam's phase front at the locations of the mirrors. 
We then use a gradient-flow minimization method to nu- 
merically search for the beam profile that corresponds 
to the minimum thermal noise, subject to the constraint 
of constant diffraction loss. While diffractional loss for 
any beam/mirror configuration will be determined from 
the clipping approximation, thermal noise will be calcu- 
lated from so-called scaling laws, which state that ther- 
mal noise spectral density is proportional to a simple 
spatial-frcqucncy-domain integral over the Fourier repre- 
sentation of the beam's intensity distribution profile on 
the mirror surface, weighted by a power of the radial wave 
number. These scaling laws apply when the mirrors are 
approximated as half- infinite (i.e., filling half of the en- 
tire space) . Their existence has been proposed by various 
authors (O'Shaughnessy, Strigin and Vyatchanin [§] and 
Vinet [§[ for substrate thermoelastic noise, Vyatchanin 
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FIG. 1: Advanced LIGO noise budget for the baseline design 
compared a hypothetical situation when the conical mirror is 
used and it only impacts coating thermal noise. This situation 
is not real because some other noises do change for better or 
for worse when one switches to the conical beam. In order to 
get the true noise curve of a LIGO-type instrument using 
conical mirrors, one would have to do a much more detailed 
analysis, well beyond the scope of this paper. 

[Io| for coating Brownian noise, and O'Shaughnessy 
for all four types of noises), while Lovelace [13] clearly 
presented and verified the scaling laws for all four types 
of noises. 

The rest of paper is organized as follows. Section II re- 
minds the reader of the Fluctuation-Dissipation-Theorem 
(FDT) -based approach to internal thermal noise due to 
Levin, and summarizes its application to the characteri- 
zation of spectral density of coating and substrate noises. 
Section III outlines the minimization problem. Section 
IV describes the results of this minimization - we discuss 
the mirror profile and the light beam that correspond to 
the minimum of substrate and coating noise. Section V 
expands on this discussion and addresses the issues of 
tolerance to imperfections and compatibility with LIGO. 
Conclusions are drawn in in Section VI. 

While we were working on this paper and after most 
results were published in the thesis of MB [l3[ , it came to 
our attention that a theoretical lower limit for the coating 
noise in an Advanced LIGO detector was derived in [l4| , 
but it was not proven that it can be reached. In this 
paper we construct an actual beam that approximately 
achieves this limit. 



I. NOISE CHARACTERIZATION 

This section briefly reviews the main results of Levin 
and Lovelace, which will be the basis of our work. 
Lovelace assumes half-infinite mirrors, i.e. he neglects 
all effects arising from a mirror's finite thickness as well 



as mirror edge effects, and also ignores the dynamics of 
the mirror, e.g., by using the quasi-static approximation. 

LIGO extracts the gravitational wave signal by mea- 
suring the position of the mirrors. The position informa- 
tion is read as q(t), a weighted average of the mirror's 
longitudinal position, which depends on Z(r, 0, t), as fol- 
lows: 

q(t)= / d<t> I drrp(r)Z(r,<f>,t) (1) 
Jo Jo 

where Z is the displacement of the mass element at (r, <f) 
of the mirror surface, R is the mirror radius, and p(r) is 
the light intensity of the axisymmetric beam at distance 
r from the optical axis, which satisfies the normalization 
of 

/ p(r)rdr = 1 . (2) 
Jo 

Internal thermal noise will cause small fluctuations in 
the longitudinal position of the mirror surface, Z(r, 0, t). 
This noise can be divided into two different types: Brow- 
nian and Thermoelastic. Brownian thermal noise is due 
to imperfections in the substrate or coating material that 
couple normal modes of vibration to each other. Ther- 
moelastic noise is due to random heat flow in the mirror 
that causes some regions to expand and some to con- 
tract. Both noises arise from the substrate as well as 
from the mirror coating. Thus, we have to deal with four 
types of noise: Coating Brownian noise, Coating Ther- 
moelastic noise, Substrate Brownian noise and Substrate 
Thermoelastic noise. 

The spectral density, S*,of the fluctuations in the mea- 
sured mirror position, q, is derived from the Fluctuation 
Dissipation Theorem using Levin's thought experiment 

_ 2fc B rw diS5 

where ks is Boltzmann's constant, T is the mirror tem- 
perature, and Wdiss is the dissipated power if a longitudi- 
nal force F is applied to the mirror surface with frequency 
/ and pressure profile p(r) , identical to light intensity pro- 
file. In initial LIGO, to keep diffraction losses under 1 
ppm per bounce, the beam radius over which 95% of the 
signal is collected is kept significantly smaller than the 
mirror radius, R, and mirror thickness, H. Lovelace was 
forced to use infinite test-mass approximation because 
the non-infinite case is too difficult to solve analytically 
and will not give a simple scaling law. He later showed 
that the infinite test mass approximation holds reason- 
ably well for beams considerably larger than in initial 
LIGO. Fused silica is found to be significantly less sus- 
ceptible to the mirror edge effects and finite effects due to 
the finite thickness of the mirror than sapphire substrate. 
By using his results, we make the same assumptions in 
computing the noise. As we will see below, the condi- 
tion that 95% of the signal is collected from a mirror 
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area smaller than R is very true in the case of Gaussian 
beam, almost true in the case of Mesa beams and almost 
false in the case of the conical beams proposed here. 

Now, because the resonant frequencies of the mirror 
are of order 10 5 Hz, far higher than the 40 — 200 range in 
question , the hypothetical force, F, can be idealized as 
quasi-static when computing the resulting strain of the 
mirror. Advanced LIGO will measure from 10 Hz to 10 
kHz; however, thermal noise is dominant only from 40Hz 
to 200 Hz. 

Thus, to compute the noise S, Lovelace [l2| substi- 
tutes in Eq. [3] the Brownian and Thcrmoclastic dissi- 
pated power, Wdiss, due to a mirror deformation with 
the same pressure distribution as p(r), the light inten- 
sity. This entire procedure is based on Levin's thought 
experiment [Til , fioj ] . This results in the following rela- 
tionships (Eq. (3.1) of Lovelace (l2|): 
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dk k n \p(k)f 



(4) 



where n = 1 for Coating Brownian and Coating Thcr- 
moclastic noise, n = for Substrate Brownian noise, 
n = 2 for Substrate Thermoelastic noise, A is a constant 
that depends on the noise type and instrumental setup 
but does not change with the beam shape, and p(k) is 
the two-dimensional Fourier transform of the power dis- 
tribution over the mirror surface, p(r): 



p(k) 



dr rJo(kr)p(r) 



(5) 



p(r) = I dk U (kr)p(k) (6) 
Jo 

In the above, Jo(x) is the 0*' 1 order Bessel function of 
the first kind. This makes numerical evaluation of these 
noises easiest in the context of a minimization code that 
requires noise to be computed a large number of different 
power profiles. In the case of coating Brownian and Ther- 
moelastic noises, we have n = 1, which converts directly 
to 
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dr r\p(r)\' 



(7) 



according to the Parseval Theorem. 

Before moving on, we mention that the thermally in- 
duced gravitational wave strain noise power, Sh{f), is 
related to S by 
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because the interferometer measures 

,(<?i- q-i) - (<?s - 94) 



(8) 



(9) 



where qi is the measured position of the i th mirror and 
L = 4 km, the arm cavity length. 



II. THE MINIMIZATION PROBLEM 

Within the paraxial approximation, Gauss-Laguerre 
beams provide a complete set of orthonormal basis vec- 
tors in the space of all possible LIGO beams propagating 
along the ±z direction, supported by cavities with axi- 
symmetry around the z axis. If we choose to center the 
cavity at z = 0, and place mirrors at z = ±L/2, then for 
infinite cavities, because of time-reversal symmetry, the 
cigenmodes will have real- valued amplitude on the z = 
plane, i.e., U(r 7 z) e R. In this way, all possible U can 
be expanded as 



U(r,z) 



OO 

E 



A n ^ n (r, z), 



(10) 



with A n <S R, and ^f n (r,z) the n-th Gauss-Laguerre 
mode with waist located at z = (and angular quan- 
tum number equal to 0): 



^n(r,z) 



(2n+l)i arctan ■ 



V2r/ 



w 



(11) 



where Wg is the waist size of the Gauss Laguerre mode, 
z = ttw'q/X, V'n( r ) = Ce~ r / 2 C n {r 2 ) and C n the n-th 
Laguerre Polynomial. The Gauss-Laguerre modes satisfy 
the following normalization: 
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^j(r, z)* fc (r, z)rdr = 6 kj , Vz . (12) 



Note that there is a freedom in choosing Gauss-Laguerre 
modes with different waist sizes, but we choose the Mini- 
mal Gauss-Laguerre mode (which has the minimum waist 
size at the mirror location), with zq = L/2. 

The infinite mirrors that support this as the fundamen- 
tal mode are determined by the constant-phase surface 
of U(r,z) around the mirror locations (z = ±L/2). The 
light intensity on the mirror surface is simply given by 



p(r) oc \U(r,z = L/2)\' 



(13) 



In the clipping approximation, the mirrors are simply 
taken as finite portions of the phase front, with r < R, 
and the corresponding diffraction loss is 

Sr 'P(r)rdr 



Jo P(r)rdr 

E,"=o ^Mj Sr Mr, zo^jr, z )rdr 



(14) 



where R = 17cm is the mirror radius in Advanced LIGO. 

Based on this, and the aforementioned scaling law for 
Si, the coating thermal noise ((?]) can thus be written as 
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a quartic function of A m : 

5ioc J2 AiAjAkAi / *i**Mj*rdr, (15) 

i,j,k,l=0 

while the normalization of p(r) requires [assuming small 
diffraction loss, Cf. Eq. ©] 

OO 

E^ = 1 ( 16 ) 

Fixing an optical loss of e per bounce adds another con- 
straint, 

°o poo 

V AiAj / *i(r, z )*i(r, z )rcfr = e . (17) 

Here we choose to have e = 10~ 6 . 

The minimization is carried out on constraint- 
satisfying sub-manifold of the space of linear combina- 
tions of Gauss-Laguerre basis functions. Thus, the gra- 
dient generally points out of the submanifold. To correct 
for this, at every step of the minimization after moving 
along the gradient of the coating noise (fll)|) , we return to 
the constraint-satisfying manifold by moving along the 
gradient of the diffraction loss untill (fl7|) is satisfied. The 
last step is renormalizing so that (|16|) is maintained as 
well. 

Numerically, we found that when minimizing (|15p sub- 
ject to the constraints ([To]) and (TTT)) in a space with a 
large number of dimensions, one runs into local minima. 
A good way [l8[ to avoid them was to strat in a space 
with few dimensions (low number of Gauss-Laguerre co- 
efficients) and increase the dimensionality one by one, 
always using the result of the previous step as innitial 
guess. 

III. RESULTS 
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FIG. 2: The three types of thermal noise for the Gaussian, 
Conical and Mesa beams normalized so that the conical beam 
noise is 1 in each cathegory 



TABLE I: Ratio of Mesa and baseline Gaussian cavity noise 
to conical cavity noise for different types of noises. 



Noise 


Mesa 


Gaussian 


Coating 


2.32 


5.45 


Substrate Thermoelastic 


3.32 


11.38 


Substrate Brownian 


1.53 


2.33 



The minimization code discussed in the previous sec- 
tion converges to a beam much wider than Mesa shown 
in Fig. [21 while maintaining the same diffraction loss 
and total power. Switching to a wider beam naturally 
leads to an overall decrease in all types of thermal noise, 
even though only the coating noise is actively minimized. 
This is illustrated in Table U and Fig. [21 The mirror pro- 
file that should be used to support this beam is shown 
in Fig. 21 Since the mirrors are approximately conical 
in shape, we will name our beam the Conical beam. In- 
terestingly, we note that the mean slope of the cone is 
roughly R/L f» 4.3 x 10~ 5 f» 0.4A/cm. Wc also note 
that the mild deviations from a perfect c one oscillates 
spatially at the Fresnel scale of y/ XL/ (27r) w 2.6 cm. 

In Fig. [31 when we compare the intensity profile of the 
Conical beam with those of Mesa and Gaussian beams, 
we found that the intensity of the Conical beam extends 



more toward boundaries of the mirrors than Mesa and 
Gaussian beams - which increases the level of averaging, 
and is critical in achieving a significantly lower thermal 
noise. This extension does not bring a heightened optical 
loss, because the Conical beam cuts off more sharply at 
the edge of the mirror than Gaussian and Mesa beams, 
as illustrated in Fig. Instead of a Gaussian-like smooth 
fall-off, power cutoff of the Conical beam is characterized 
by oscillations outside the mirror radius. We attribute 
such an oscillatory cutoff to the mild oscillations of mirror 
surface shown in Fig. 01 Overall, because of extending 
to larger radii, given the same optical power, the light 
intensity of a Conical beam everywhere is less then the 
central area of the baseline Gaussian, while only in the 
central peak does the Conical beam slightly surpasses 
the Mesa plateau value. As a consequence, our beam 
should be comparable to Mesa in terms of field-strength 
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FIG. 3: The power distribution of the lowest-noise 35 coef- 
ficient beam found using our minimization algorithm, com- 
pared with the previously published Mesa Beam. 




FIG. 4: Mirror height, plotted in units of A = 1.06/xm, for 
the nearly conical (solid blue line) mirror compared with the 
its nearly concentric (dashed red line) and nearly flat (dotted 
green line) Mesa counterparts. 



tolerability by coating materials. 

It is worth mentioning that conical cavities similar to 
ours have been proposed for the generation of Bessel- 
Gauss beams pj|. However, our beam cuts off more 
sharply than Bessel-Gauss beams at the edge of the mir- 
rors. This must also be due to our mirrors' oscillatory 
deviations from perfect cones. 



A. Convergence and Conical Cavities with Fewer 
Coefficients 

The results presented throughout most of this paper 
are based on our 35 coefficient minimization result. Low- 
noise cavities may however be found that employ fewer 
Gauss-Laguerre coefficients. They are not as good as our 
final result, but still much better than Mesa. Employing 
fewer Gauss-Lagucrrc modes may lead to a cavity that 
is more desirable for reasons other than thermal noise. 
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FIG. 5: The power distribution of the 35 coefficient lowest- 
noise beam outside the mirror compared to the theoretical 
prediction for Mesa. In the clipping approximation, the inte- 
gral of this power is the assumed to be the diffraction loss. 
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FIG. 6: Power distribution at the mirror as a function of the 
number of Gauss-Laguerre coefficients employed in the mini- 
mization code. When fewer coefficients are used, the beam is 
spread over a smaller area of the mirror. 



Another reason for showing this figure is to demonstrate 
convergence. 

When using fewer Gauss-Laguerre modes in the beam 
expansion, the beam does not extend all the way out to 
the end of the mirror. In Fig. [51 we show the intensity 
profile at the mirror for beams obtained by minimizing 
the coating noise over the lowest 7, 14, 21, 27 and 35 
Gauss-Laguerre modes. As seen in Fig [5J there is little 
qualitative difference in the mirror shapes, depending on 
the number of coefficients used. 
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FIG. 7: The difference in the power distribution at the mirror 
between the 35 coefficient mode and several other modes with 
a different number of Gauss-Laguerre coefficients employed in 
the minimization code. 



FIG. 9: Absolute Value of eigenvalues of the one-way propa- 
gator. The blue dots represent n = or axisymmetric modes, 
the green n = 1, and the red n — 2. 



mirror to the other: 
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FIG. 8: Mirrors supporting the 7, 12, and 30 coefficients 
modes. The mirror height is measured in units of A, where 
A = 1.06/xm is the wavelength of the light used in the inter- 
ferometer. 



(18) 

where J m (z) is the m th order Bcssed function of the first 
kind given by 



Jn(z) 



1 



2?ri r 



iz cos (p in<fi 



dcj) 



(19) 



L is the length of the arm cavity, k is the wave num- 
ber k — and hi^{r) are the mirror heights and are 
assumed to be equal. For the conical mirror shape ob- 
tained from the previous section, we then computed the 
eigenvalues of the axisymmetric propagator and obtained 
only one that was close to 1. All other higher eigenmodes 
had high diffraction losses. Figure [9] we show the ab- 
solute values of the eigenvalues, for modes with n — 
(axisymmetric), 1 and 2. The mode with minimum loss 
has |Ai| = 1 — 1.45 • 10~ 6 corresponding to a per-bounce 
diffraction loss of 



1 - |Af | = 2.9ppm. 



(20) 



B. Optical Modes Supported by Finite 
Nearly-Conical Mirrors and their Diffraction Loss 



We designed our mirrors to have diffraction losses of 1 
ppm in the clipping approximation in order to agree with 
Mesa and the Baseline Gaussian designs previously con- 
sidered [2(j ■ Here we need to verify that the fundamental 
mode supported by a finite conical cavity with radius R is 
indeed close to the beam we constructed, with a diffrac- 
tion loss close to 1 ppm. We are also interested in general 
the diffraction loss of higher modes. 

Following [2lT ] , for modes with angular quantum num- 
ber n, we construct a 1-D radial propagator from one 



This is much more than the 1 ppm we find in the clip- 
ping approximation and illustrates the limitations of that 
method. However, this does not invalidate our configura- 
tion because 3 ppm is still reasonable for use in Advanced 
LIGO, and when the per-bounce losses of Gaussian and 
Mesa beams are relaxed to 3 ppm, their thermal noises 
only change by a very small amount. 



IV. TOLERANCE TO IMPERFECTIONS AND 
COMPATIBILITY WITH LIGO 

To investigate the effects of mirror perturbations on 
the conical cavity, we simulated the 2-D propagation of 
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FIG. 10: Diffraction loss in a conical cavity as a function of 
the tilt, when its mirrors are symmetrically tilted. 



the light from one mirror of the cavity to the other, using 
an FFT code [22j|, without assuming any symmetry. We 
used this tool to study the effects of mirror tilt, mirror 
displacement, and mirror figure errors as well as checking 
the diffraction losses of the conical cavity. During light 
propagation in the cavity, we note that reflection from 
the mirror is a diagonal operator in position space, while 
propagation in free space from one plane to the other 
is diagonal in spatial- frequency domain. Therefore the 
most efficient means of propagation is as follows. When 
propagating between planes at the mirror locations, use 
spatial-frequency domain, and when reflecting from mir- 
ror surfaces, use position domain — and insert 2-D FFT 
between these processes (and hence the name FFT code) . 
For a given cavity, we propagate and initial field multiple 
times in order to obtain the mode with the lowest loss. 
As a test, we first used the code on the calculated conical 
mirror shape, and obtained a beam very close to the theo- 
retical prediction. The diffraction loss per bounce is 3.03 
ppm, which agrees well with (|20[) . In the following, we 
use the FFT code in various situations of imperfection, 
namely mirror tilt, translation and mirror figure error, 
and estimate the requirement on these imperfections, if 
the conical cavity is to be used in LIGO. 



A. Sensitivity to Mirror Tilt 

Our numerical simulation shows Mesa cavities to in- 
crease the diffraction loss to about 3 ppm when perturbed 
symmetrically by a 10~ 8 radians mirror tilt. The same 
tilt induces a diffraction loss of about 70 ppm in the con- 
ical cavity. (Antisymmetric tilt can be treated as trans- 
lation of one mirror, which will be considered next.) The 
diffraction loss depends quadratically on the tilt angle as 
shown in Figure fTUl The interpolating function, when x 



FIG. 11: Dots represent the diffraction loss in ppm for a con- 
ical cavity perturbed by moving one mirror away from the 
optic axis as a function of this displacement. The continuous 
line is a quadratic function fit to the data. 

is the tilt in radians, is given by 

loss = ( 3.03 + 0.646^—^) ppm (21) 
\ nrad^ / 

As a consequence in order for diffraction loss not to in- 
crease significantly due to tile error, the conical mirrors 
proposed here respond to tilt perturbation more strongly 
than Mesa does. The suspension and mirror control sys- 
tem would need to be engineered to control the mirror 
direction better. If LIGO is to use conical mirrors, the 
tilt needs to be controlled at the level of about 3 nano 
radians, in order for diffraction loss due to tilt error not 
to exceed 10 ppm. 

B. Sensitivity to Mirror Translation 



Simulations show that a conical cavity with a mirror 
translated by 10 /xm from its ideal position will have a 
diffraction loss of 57.61 ppm. A similar diffraction loss 
is seen in a Mesa cavity with a 4 mm error in the mirror 
positioning. Thus the conical cavity is far more sensitive 
to mirror translation than Mesa. A quadratic function 
can be fit well to the as seen in Fig. [TlJ The diffraction 
loss in parts per bounce per million of a conical cavity 
perturbed by moving one of the mirrors a distance x par- 
allel to itself in a direction orthogonal to the optic axis 
is thus approximativcly given by the formula 

loss = ( 3.03 + 0.542 ppm (22) 
y /inr J 

As a consequence, if Advanced LIGO is to use conical 
mirrors, mirror position needs to be controlled to about 
4 /urn, in order for diffraction loss due to displacement 
error not to exceed 10 ppm. If one uses antisymmetric 
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mirror tilt to compensate errors in mirror location, the 
tilt would need to be controlled to about 4 /jm/4km = 1 
nrad. 



C. Mirror Figure Error and Contribution from 
Different Scales 

The mirror figure errors are deviations of the mirror 
surface heights from their theoretical values. Figure er- 
rors have been measured experimentally for the LIGO I 
mirrors currently used in the experiment. We used the 
real measured data from [23|. Since LIGO I and Ad- 
vanced LIGO have different mirror sizes, we interpolated 
the data and stretched it from 12.5 cm (LIGO-I) to 17 cm 
(Advanced LIGO), in the same way as Ref. [1JJ. (This 
will not be very realistic, because the length scale of the 
perturbations will make a difference in the losses.) Af- 
ter numerically solving for the lowest-loss mode using the 
FFT code, we found that Conical cavities are much more 
sensitive to this type of figure error, by giving a loss of 
405 ppm, than Mesa cavities which gives a loss of 5 ppm. 
Nevertheless, if we demand a figure error of 1/10 that of 
the LIGO-I stretched error, we would recover a loss of 
5 ppm, which then becomes reasonable. 

In order to guide further development of mirror man- 
ufacturing, it is interesting to study the contribution 
to optical loss from figure errors at different scales. 
Large-scale, or low-spatial-frequency errors cause light to 
slightly deviate from the cavity axis, while small-scale, or 
high-spatial-frcqucncy errors will cause light to deviate 
more significantly from the axis. In this way, for both 
Mesa and conical cavities, errors at high spatial frequen- 
cies will generate loss anyway, while for Mesa cavities, 
low-spatial-frequency errors can be less dangerous, since 
Mesa beam does not cut off so sharply at the edges of 
the mirrors. The other way to look at it is that unlike 
the Conical cavity, the Mesa cavity has more than one 
low-diffraction-loss modes, which provides a "reserve" 
to maintain a low-loss fundamental mode. The conical 
cavity's sensitivity to tilt and displacement agrees with 
this argument, because tilt and displacement can be con- 
sidered as low- frequency imperfections. Following this 
argument, we also note that since the LIGO-I error is 
stretched, it tends to decrease the spatial frequency of 
the figure error, therefore will be biased against conical 
cavities. 

Now we numerically study contributions from differ- 
ent spatial scales by applying high-pass spatial filters to 
the stretched LIGO-I noise, and then compute diffrac- 
tion loss of the resulting fundamental optical mode. As 
we see from Table HH as we gradually allow low-spatial- 
frequency errors to enter, the Conical and Mesa losses 
initially trace each other — until a particular spatial 
scale (in our case A ~ R), when Conical loss increases 
dramatically, eventually climbing up to 100 times that of 
the Mesa loss, while Mesa loss first increases only mildly, 
and then even decreases, probably due to particular char- 



maximum A 
included 


Cone 


Mesa 


R/16 


3.16 


2.68 


R/8 


4.71 


3.83 


R/4 


6.34 


5.24 


R/2 


12.47 


7.96 


R 


56.28 


5.63 


2R 


184.13 


5.03 


4R 


404.83 


4.86 



TABLE II: Optical losses in ppm of conical and mesa cavi- 
ties, with LIGO-I mirror figure error stretched and high-pass 
filtered. When bandwidth of the filter is increased, including 
more low-frequency (i.e., long length-scale) fluctuations, the 
loss of the conical cavity increases monotonically, while that 
of the mesa cavity first increase and then slightly decrease. 

acteristics of this data set. This means that, in order to 
make a conical mirror within the LIGO specification, we 
only need to focus on large-scale errors (A > R) — while 
roughness at smaller scales at LIGO-I level is already ac- 
ceptable for conical cavities. 



V. CONCLUSIONS AND FURTHER 
DISCUSSIONS 

We developed and implemented a simple minimization 
algorithm and used it for finding the LIGO-compatible 
beam with the lowest thermal noise. The result of this 
minimization is a beam with thermal noise a factor of 
2.32 (in power) lower than previously considered Mesa 
Beams and a factor of 5.45 (in power) lower than the 
Gaussian beams employed in the current baseline design. 
The mirror that supports this mode is found to have 
nearly conical shape. Using an axisymmetric 1-D prop- 
agator, we found that contrary to spherical and Mesa 
cavities, the conical cavity only has one cigenmode with 
very low diffraction loss (the fundamental mode), while 
higher modes have much higher optical losses. 

By using an FFT propagation code, we have analyzed 
the conical cavity's tolerance to practical imperfections. 
Qualitatively, the conical cavity is much more suscepti- 
ble than Mesa cavities, to imperfections with low spatial 
frequencies — including tilt and mirror translation. At 
higher spatial frequencies, the conical cavity is compara- 
ble to Mesa cavities. This behavior can be explained by 
the fact that the presence of additional low-loss modes 
allows the Mesa cavity to continue having low diffrac- 
tion loss when lower-spatial-frcqucncy perturbations arc 
made to the mirrors. 

Recently, Parametric Instability was shown to be a se- 
rious problem in Advanced Gravitational Wave Detectors 
[24l I25I l26l . [27T | . It arises when the beat frequency of two 
optical modes is close to the mechanical frequency of an 
acoustic mode of the mirror. Wc believe Conical Beams 
exhibit lower Parametric Instability then Mesa because 
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our higher order modes are very lossy and little power 
is available to excite the mechanical modes. The topic 
needs to be thoroughly researched in the future. 

Our numerical simulations show that, in order to 
achieve diffraction losses close to Mesa on a Conical Cav- 
ity, the mirror needs to be manufactured or corrected 
with a CO2 Laser such that its large scale deviations 
from the desired shape are roughly 10 times lower then 
in the case of Initial LIGO. The mirror orientation needs 
to be controlled to about 1 nanoradian which corresponds 
to controlling the mirror displacement to the accuracy of 
several microns. 

Even though the Conical attains the lowest possible 
coating thermal noise, compromise configurations like 
those considered in [28[ will be found to provide a rea- 
sonable noise level while being easier to implement. 

Acknowledgments 

We thank Kip Thorne for suggesting this problem to 
us. We thank GariLynn Billingsley for providing infor- 
mation about mirror manufacturing and Rana Adhikari 



for information about suspension systems of future LIGO 
detectors. We thank Jayashree Balakrishna, Gregory 
Daues, Ruxandra Bondarcscu, Andrew Lundgren and 
Dave Tsang for carefully proofreading the manuscript, 
checking our calculations, and valuable advice on the 
content and the format of this paper. We thank Kentaro 
Somiya for helping us with Fig. [TJ M.B. and Y.C. were 
supported by the Alexander von Humboldt Founda- 
tion's Sofja Kovalevskaja Programme, NSF grants PHY- 
0653653 and PHY-0601459, as well as the David and Bar- 
bara Groce startup fund at Caltech. MB is grateful for 
travel support generously offered by The Center for Com- 
putation and Technology at The Louisiana State Univer- 
sity, Rochester Institute of Technology, The American 
Physical Society, The University of Tuebingen, Univer- 
sita Degli Studi Del Sanio, and Jayashree Balakrishna. 
While preparing this paper, we enjoyed inspiring and 
stimulating conversations with Gabrielle Allen, David 
Blair, Manuela Campanelli, Erika D Ambrosio, Vincenzo 
Galdi, Kazuaki Kuroda, Konstantinos Kokkotas, Carlos 
Lousto, Vincenzo Pierro, Innocento Pinto, Ed Seidel, Du- 
mitru Vulcanov and Hiro Yamanoto. 



[1] 
[2] 
[3] 
[4] 
[5] 



[6] 



[9] 



[10] 



[11] 

[12] 
[13] 



http://www.ligo.caltech.edu/! 
http:/ /www.virgo. infh.it/ 1 
http:/ /www. geo600.de/ 
http://tamago.mtk.nao.ac.jp/ 

G. M. Harry, H. Armandula, E. Black, D. R. M. Crooks, 

G. Cagnoli, J. Hough, P. Murray, S. Reid, S. Rowan, P. 
Sneddon, M. M. Fejer, R. Route, and S. D. Penn, Appl. 
Opt. 45, 1569-1574 (2006) 

H. Mueller-Ebhardt, H. Rehbein, R. Schnabel, K. Danz- 
mann, Y. Chen, Phys. Rev. Lett. 100, 013601 (2008). 

T. Corbitt, Y. Chen, H. Mueller-Ebhardt, E. Innerhofer, 
D. Ottaway, H. Rehbein, D. Sigg, S. Whitcomb, C. Wipf, 
N. Mavalvala Phys. Rev. Lett. 98, 150802 (2007). 
J. Vinet, Class. Quantum Grav. 22, 1395 (2005). 
R. O'Shaughnessy, S. Strigin, S. Vyatchanin, "The 
implications of Mexican-hat mirrors: calculations 
of thermoelastic noise and interferometer sensitivity 
to perturbation for the Mexican-hat-mirror proposal 
for advanced LIGO," (submitted to Phys. Rev. D), 
arXiv:gr-qc/0409050| 

Vyatchanin. "Fdt approach calcula- 
tions of brownian noise in thin layer" 
(internal LIGO document t040242-0q z) , 
|http: / /www.ligo.caltech.edu /docs /T /T040242-00.pdf | 
R. O'Schaughnessy, Class Quant.Grav. 23 (2006) 7627, 
|arXiv:gr-qc/0607035 



G. Lovelace, Class. Quant. Grav. 24, 4491 (2007). 
Mihai Bondarescu, "Topics in General Rel- 
ativity", Ph.D. Thesis, California Institute 
of Technology, Pasadena, California, 2007. 



|http://etd.caltech.edu/etd/available/etd-05282007-231321/ 
[14] V. Galdi, G. Castaldi, V. Pierro, I. M. Pinto, J. Agresti, 
E. D'Ambrosio and R. DeSalvo, Phys. Rev. D 73, 127101 
(2006). 



[15 
[16 



[17 
[18 
[19 
[20 
[21 

[22 

[23 
[24 

[25 
[26 

[27 
[28 
[29 



Jr and J.H. Eberly, Phys. Rev. 
S. Thorne, Phys. Rev. D 74, 



Y. Levin, Phys. Rev. D 57, 659 (1997). 
Kip S Thorne, Mihai Bondarescu, Yanbei Chen "Grav- 
itational Waves - A Web Based Course," Available on- 
line at |http://elmer.tapir.caltech.edu/ph237| or on DVD 
from Caltech library. The Levin Thought Experiment is 
discussed in the week 12 homework. 
M. Born and E. Wolf, Principles of Optics, 2nd edition, 
New York, The MacMillan Company, 1964. 
M. Bondarescu, M. Alcubierre and E. Seidel, Class. 
Quant. Grav. 19, 375 (2002). 
J. Drunin, J. J. Miceli, 
Lett. 58, 1499 (1987). 
M. Bondarescu and K. 
082003 (2006). 

P. Savov, J. Agresti, Y. Chen, E. D'Ambrosio, "A 
duality relation between non-spherical optical cavities 
and its application to gravitational-wave detectors," 
arXiv.org jgr-qc/0511062) 

see, e.g., Erika DAmbrosio, Phys. Rev. D 67, 102004 

(2003) 

LIGO http: //www.ligo.caltech.edu/e2e l 
V.B. Braginsky, S.E. Strigin and S.P. Vyatchanin, 
Phys. Lett. A 287, 331 (2001); V.B. Braginsky, S.E. Stri- 
gin and S.P. Vy atchanin, Phys. Lett. A 305, 111 (2002). 
C. Zhao et al. larXiv:0801.1150l [gr-qc], 
C. Zhao, L. Ju, H . Miao, S. Gras, Y. Fan and D. G. Blair, 
larXiv:0710.2383l [gr-qc] . 

C. Zhao, L. Ju, J. Degallaix, 

D. G. Blair, Phys. Rev. Lett. 94 
[arXiv:gr-qc/0502079] . 
A. P. Lundgren, R. Bondarescu, 
M. Bondarescu, Phys. Rev. D 77, 



S. Gras and 
121102 (2005) 



D. Tsang and 
042003 (2008) 
|arXiv:0710.3808l [gr-qc]] . 

for example, a decrease in total noise by a factor of two 



10 



corresponds to an increase in event rate of observing grav- 
itational waves roughly by a factor of three 



